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Abstract 

A mean field spin system consisting two interacting groups each with homogeneous 
interaction coefficients is introduced and studied. Existence of the thermodynamic 
limit is shown by an asymptotic sub-addittivity method and factorization of correla- 
tion functions is proved almost everywhere. The free energy solution of the model is 
obtained by upper and lower bounds and by showing that their difference vanishes 
for large volumes. 

1 Introduction 

In this work we consider the problem of characterizing the equilibrium statistical mechan- 
ics of an interacting system of two set of spins. We aim to tackle the most general case 
of a two population system in the mean field approximation. 

Mean field two-population models have been useful since the study of metamagnets 
which started in the '70s (see [HIS]), and have been encountered very recently in the study 
of loss of gibbsianness for a system whose evolution is described by Glauber dynamics (see 
[3]). In Refs [H E] a two-population mean-field model is used as an approximation to a 
bipartite lattice assumed to describe an antiferromagnetic system, and is found to repro- 
duce qualitatively the expected phase transitions, which are then studied at criticality. 
In Ref. [3], instead, particles are subject to a time-evolving random field which acts on 
particles by partioning them in two groups, leading to a mean field model mathematically 
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analogous to the former. In this case a characterization of the whole phase diagram is 
provided. 

The systems considered by both works can be seen as a restriction of a more general 
model, which is presented here, and which arises naturally as an interacting mixture of 
two systems of Curie- Weiss type. 

Our results can be summarised as follows. After introducing the model we show in 
section 3 that it is well posed by showing that its thermodynamic limit exists. The result 
is non trivial because sub-additivity is not met at finite volume. In section 4 we show 
that the system fulfills a factorization property for the correlation functions which reduces 
the equilibrium state to only two degrees of freedom the equilibrium state. The method 
is conceptually similar to the one developed by Guerra in [4j to derive identities for the 
overlap distributions in the Sherrington and Kirkpatrick model. 

We also derive the pressure of the model by rigorous methods developed in the recent 
study of mean field spin glasses (see [6] for a review). It is interesting to notice that though 
very simple, our model encompasses a range of regimes that do not admit solution by the 
elegant interpolation method used in the celebrated existence result of the Sherrington and 
Kirkpatrick model [7j. This is due to the lack of positivity of the quadratic form describing 
the considered interaction. Nevertheless we are able to solve the model exactly, section 
5, using the lower bound provided by the Gibbs variational principle, and thanks to a 
further bound given by a partitioning of the configuration space, itself originally devised 
in the study of spin glasses (see P, El E])- 

As in the classical Curie- Weiss model, the exact solution is provided in an implicit 
form; for our system, however, we find two equations of state, which are coupled as well 
as trascendental, and this makes the full characterization of all the possible regimes highly 
non-trivial: a robust numerical analysis becomes essential and can be found in [9j, where 
an application to social sciences is considered. 

Some aspects of the regimes can nonetheless be studied analytically, and this is done 
in section [6} while a global study of the phase diagram for our model is left to be carried 
on in a future work. 
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2 The Model 



Our model is defined by the Hamiltonian 



H 



1 ^ 



(2.1) 



We consider Ising spins, (Tj = ±1, and symmetric interactions Jjj. We divide the 
particles P = {1, 2, 3, A^} into 2 types A and 5 with A U E = P, A n S = 0, and sizes 
iVi = |A| and A^2 = \B\, where A'^i + A^2 = N. Given two particles i and j, their mutual 
interaction parameter J^j depends on the subset they belong to, as specified by the matrix 

Ni N2 
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where each matrix block has constant elements: J\\ and J22 tunc the interactions within 
each of the two types, and J12 controls the interaction between two particles of different 
types. In view of the applications considered in the introduction, we assume Jn > and 
J22 > 0, whereas J12 can be either positive or negative. 

Analogously, the field hi takes two values h\ and h^^ depending on the type of i, as 
described by the following vector: 

A^i I / hi \ 



iV2 < 

By introducing the magnetization of a subset S as 



ms{(y) 



-T 



1^1 



i<=S 



and indicating by mi and m2 the magnetizations within the subsets A and B and by 
— the relative size of subset A on the whole, we may easily express the Hamiltonian 



per particle as 

^ ^ — —- Ji2a(l — Q;)mim2 + J22(l — tt)^^i2] —h\ami — h2{l — a)m2 (2.2) 

The usual statistical mechanics framework defines the equilibrium value of an observ- 
able /(cr) as the average with respect to the Gihhs distribution defined by the Hamiltonian. 
We call this average the Gihhs state for /(a), and write it explicitly as: 

The main observable for our model is the average of a spin configuration, i.e. the 
magnetization, m{a), which explicitly reads: 

N 

N 



i=l 



Our quantity of interest is therefore (m): to find it, as well as the moments of many 
other observables, statistical mechanics leads us to consider the pressure function: 



N 

a 

It is easy to verify that, once it's been derived exactly, the pressure is capable of 
generating the Gibbs state for the magnetization as 

In order to simplify the analytical study of the model it is useful to observe that the 
Hamiltonian is left invariant by the action of a group of transformations, so that only a 
subspace of parameter space needs to be considered. 

The symmetry group is described by G = Z2 x Z2 x Z2. 

We can represent a point in our parameter space as (m, J, h, a), where 

mi 

m=| |,n=| \ , a 






1712 

Therefore, given the limitations on the values of our parameters, the whole parameter 
space is given hy S ^ [-1, 1]^ x x x x [0, 1]. 



If we consider the representation of G given by the 8 matrices 




+1 or — 1 and 




?7j = +1 or — 1 



we can immediately reahze that G is a symmetry from the Hamiltonian representation: 



3 Existence of the thermodynamic hmit 

We shall prove that our model admits a thermodynamic limit by exploiting an existence 
theorem provided for mean field models in p!0]: the result states that the existence of 
the pressure per particle for large volumes is guaranteed by a monotonicity condition on 
the equilibrium state of the Hamiltonian. Such a result proves to be quite useful when 
the condition of convexity introduced by the interpolation method [TJ E] doesn't apply 
due to lack of positivity of the quadratic form representing the interactions. We therefore 
prove the existence of the thermodynamic limit independently of an exact solution. Such 
a line of enquiry is pursued in view of further refinements of our model, that shall possibly 
involve random interactions of spin glass or random graph type, and that might or might 
not come with an exact expression for the pressure. 

Proposition 3.1 There exists a function p of the parameters (a, Ji,2, </2,2, ^i, ^2) 
such that 



The previous proposition is proved with a series of lemmas. Theorem 1 in pO] states that 
given an Hamiltonian Hj^ and its associated equilibrium state the model admits a 
thermodynamic limit whenever the physical condition 




lim pn = p . 





(3.3) 



is verified. 
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We proceed by first verifying this condition for a working Hamiltonian ifjv, and then 
showing that its pressure pN tends to our original pressure as N increases. We choose 



if TV in such a way that the condition (3.3) is verified as an equahty. 



Our working Hamiltonian Hn is defined as follows: 



fr _ fV(l) I f>{12) , fr{2) 



where 



i^j=l...Ari ^ i^j=l...N2 



iV 

1=1. ..JVx 

Lemma 3.1 There exists a function p such that 

lim pN = P 

Proof: By definition of H^^ and by the invariance ofooN with respect to spin permutations, 

{aN - l)aN 

We can find a similar form for if and H^^'^\ which implies that for any two positive 
integers Ki + K2 = N we have 



which verifies (3.3) and proves our Lemma. □ 



The following two Lemmas show that the difference between and if a? is thermody- 
namically negligible and as a consequence their pressures coincide in the thermodynamic 
limit. 

For convenience we shall re-express our Hamiltonian in the following way: 

TT _ tt{1) , rr(12) , rr(2) 



where we define 



i,0=\...Nx i,j=l...N2 



i=l...Ari 
j=l...JV2 



Lemma 3.2 

i^jv = ^jv + 0(l) (3.4) 



lim — — = lim — — 

N-^oo N Ar-»oo N 



Proof: 



i,j=l...Ni " ^ i^j=l...Ni " i=l...Ni 



and since — jy 



1 

..N(aN - 1) ^ 



aN-1 ^ 

T^OiJ 

aN 

i^j=l...Ni 



and so 

//« = ^« + 0(l). 

We can similarly get estimates for H^^ and H^"^^ in terms of and H^'^\ which 
implies 

□ 

Lemma 3.3 Say pn — — InZjv, and say /iAr(cr) = — ^^fi''^^ ^> ^^'^ ^'^ ^'^ 
analogous way. 
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Define 



Then 



kN = \\hN — hN\\ = sup {\hj^{a) — hj^{a)\} < CO. 
<7e{-i,+i}^ 



(3.5) 



\Pn-Pn\ < \\hN - h 



N\ 



Proof: 



Pn-Pn 



— Ill Zat 111 Zm = — In 

N ^ N ^ N Z 



N 



— In < — In 



N 



In 



-Nkf 



AT 



where the inequahty follows from the definition of k]\j in (3.5) and from monotonicity of 



the exponential and logarithmic functions. The inequality for — pn is obtained in a 
similar fashion. □ 



We are now ready to prove the main result for this section: 



Proof of Proposition 3.1: The existence of the thermodynamic limit follows from our 



Lemmas. Indeed, since by Lemma 3.1 the limit for pjv exists, Lemma 3.3 and Lemma 3.2 
tell us that 

lim \pn — Pn\ ^ lim II/iat — /iatII = 0, 



implying our result. 



□ 



4 Factorization properties 

In this section we shall prove that the correlation functions of our model factorize com- 
pletely in the thermodynamic limit, for almost every choice of parameters. This implies 
that all the thermodynamic properties of the system can be described by the magneti- 
zations nil and m2 of the two subsets A and B defined in Section |2} Indeed, the exact 
solution of the model, to be derived in the next section, comes as two coupled equations 
of state for mi and m2. 
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Proposition 4.2 



lim [uN^aiCTj) - ujNi(^i)^^N{<7j)) = 



TV— »oo 



for almost every choice of parameters {a, Jn, J12, J22, hi, , where ai, aj are spins of 
any two particles in the system. 

Proof: We recall the definition of the Hamiltonian per particle 

— = ~2 [-'^iitt^"^i+2Ji2tt(l — Q;)'T^i'Ti2 + -^22(l — tt)^"^2] ~ hiarrii — h2{l — a)m2, 
and of the pressure per particle 



Pn 



N ^ 



By taking first and second partial derivatives of p^ with respect to hi we get 

dpN 1 V- , t ^ ^'P^ 

y^aJMmi[a) — - — = aujN['mi), 



dhi Zn " ' dhf 

By using these relations we can bound above the integral with respect to hi of the 
fluctuations of mi in the Gibbs state: 



(2) 



(1) 



(cjjv(m^) — Co'Ar(mi)^) dhj 



1 



1 



h^^' 02 
\w dhl 



dhi 



1 



1 dpN 



dhl 



(1) 



Na 



{\^N{mi\(2)\ + \ujN{mi\m\) = 0{^). 



(4.6) 



On the other hand, given any a G (0, 1) we have that 

1 dpN 



ujN{mi) 



and 



a dhl ' 

2 dpN 
dJii ' 



so, by convexity of the thermodynamic pressure p = lim pj^, both quantities 

N->oo ohi 

dpN 



and 



dJii 

implies that 



have well defined thermodynamic limits almost everywhere. This together with (4.6) 
lim (co'Ar(m^) — UNimiY) = a.e. in hi, Ju. (4.7) 

JV— >oo 



In order to prove our statement we first consider spins of particles of type A, which 
we shall call ^j. Translation invariance of the Gibbs measure tells us that 
1 



i=l 

N- 



UJn\ 



1 1 1 - 1 



i=j=l 



We have that (4.8) and (4.7) imply 



lim ujNi^iQ - ^^Ni^i)uJN{^j) = 0, 



(4.9) 



which verifies our statement for all couples of spins i ^ j oi type A as defined in section 



2] (the case i = j verifies (4.9) trivially). 

Working in strict analogy as above we also get 



lim {uN{m\) — uj]y{'m2Y) = a.e. in h2, J; 



22- 



(4.10) 



Furthermore, by defining Variv(mi) = [ujj\f{ml) — co'iv(mi)^), and analogously for m2, we 



exploit (4.7) and (4.10), and use the Cauchy- Schwartz inequality to get 

\ujN{mim2) - co'Ar(mi)cjAr(m2)| ^ A/Var(mi)Var(m2) — > a.e. in Jn, J^, J22, hi, h2 

(4.11) 



By using (4.10) and (4.11) we can therefore verify statements which are analogous to 



4.9, but which concern UNi^iVj) and uj^irjirij) where ^ are spins of type A and rj are spins 



of type B. 

We have thus proved our claim for any couple of spins in the global system. 



□ 



5 Solution of the model 

We shall derive upper and lower bounds for the thermodynamic limit of the pressure. 
The lower bound is obtained through the standard entropic variational principle, while 
the upper bound is derived by a decoupling strategy. 
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5.1 Upper bound 

In order to find an upper bound for tlie pressure we sliall divide tlie configuration space 
into a partition of microstates of equal magnetization, following [HI El [5] . Since subset 
A consists of Ni spins, its magnetization can take exactly Ni + 1 values, which are the 
elements of the set 



R 



Ni 



1-1,-1 + — ,.. . ,1 —,l]- 



Clearly for every mi (a) we have that 

where 6x,y is a Kronecker delta. We can similarly define a set Rn2, so we have that 
Zn = ^exp { ^(Jiia^mi + 2Ji2a(l - a)mim2 + J22(l - a)^"^2) + /iiA/'imi + /i2A^2'"^2} 
= ^ ^ 5mi,Mi'^m2,M2exp { ^(Jiia^m^ + 2Ji2a(l - a)mim2 + J22(l - tt)^"^2) + 

+ hiNimi + h2N2m2] . (5. 

Thanks to the Kronecker delta symbols, we can substitute nii (the average of the 
spins within a configuration) with the parameter (which is not coupled to the spin 
configurations) in any convenient fashion, and the same holds for m2 and /i2- 

Therefore we can use the following relations in order to linearize all quadratic terms 
appearing in the Hamiltonian 





= 0, 


-/i2)^ 


= 0, 


- ^^2) 


= 0. 



Once we've carried out these substitutions into (5.12) we are left with a function which 
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depends only linearly on mi and 777.2: 



Zn = ^ ^ Smi,fiiSm2,n2 exp {^(JiiQ;^(2r77i//i - III) + 2Ji2Q;(1 - a)miiJ,2 + 



+ </22(l - Q;)^(27772/72 - jjl)) + 2Ji2a{l - a)7772/7l + 2Ji2Q;(l - Q;)/7i/72 + 



+ hiNimi + h2N2m2], 



and bounding above the Kronecker deltas by 1 we get 

^ 

~2 



-Z'at ^ ^ ^ exp { - y(JiiQ;Vi + 2Ji2Q;(1 - Q;)Ati//2 + ^22(1 - q;)V2) + 



+ ((JiiQ;/7i + Ji2(l - Q;)/i2 + hijNimi + (Ji2a;/7i + J22(l - Q;)/72 + /i2)iV2"^2) }• 

(5.13) 

Since both sums are taken over finitely many terms, it is possible to exchange the order 
of the two summation symbols, in order to carry out the sum over the spin configurations, 
which now factorizes, thanks to the linearity of the interaction with respect to the ms. 
This way we get: 



Zn ^ ^ G(/7i,/72). 



where 



G{iii,ii2) = exp { - y(JiiQ;Vi + 2Ji2a(l - Q;)//i//2 + ^22(1 - Q;)V2)} • 

•2^1 (cosh (Jiiq;/7i + Ji2(l - a)ii2 + /ii))'^'2^^" (cosh(Ji2Q;/7i + J22(l - a)ii2 + h2))^^ 

(5.14) 

Since the summation is taken over the ranges Rn^ and -Ratj, of cardinality A^i + 1 and 
N2 + 1, we get that the total number of terms is (A^i + 1)(A^2 + !)■ Therefore 

Zn < (iVi + l)(iV2 + l) supG, (5.15) 



which leads to the following upper bound for Pjv: 



Pn = ^InZjv^ ^ln(7Vi + l) + ^ln(iV2 + l) + ^lnsupG. (5.16) 
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Now defining the independent function 

Pc/p(/ii,/Lt2) = ^ InG = ln2 - ^(Jiia^i + 2Ji2a(l - a)/ii/U2 + J22(l - a)V2) + 

+alncosli(Jiia/ii + Ji2(l — a)/i2 + ^i)) + 

(1 - a) lncosli(Ji2Q;/xi + ^22(1 - C()n2 + ^2)), (5-17) 

tlie tliermodynamic limit gives: 

limsupP/v^ sup p(7p(/ii, /i2). (5.18) 

We can summarize the previous computation into the following: 



Lemma 5.4 Given a Hamiltonian as defined in (2.2), and defining the pressure per par- 



ticle as Pn = 4 In given parameters Jn, J12, J225 hi, /i2 and a, the following inequality 



N 

holds: 



where 



limsuppAT ^ sup Pup 



Pup = ln2 - -(Jiitt^/ii + 2Ji2a(l - a)/ii/i2 + ^22(1 - a)V2) + 
+Q;lncosh(Jiia;yUi + Ji2(l — «)yU2 + hij) + 

(1 - a) lncosh(Ji2a/ii + J22(l - «)/i2 + /i2)), (5.19) 
and (/ii,/U2) e [-1, 1]^. 

5.2 Lower bound 

The lower bound is provided by exploiting the well-known Gibbs entropic variational 
principle (see [TT], pag. 188). In our case, instead of considering the whole space of 
ansatz probability distributions considered in [H], we shall restrict to a much smaller 
one, and use the upper bound derived in the last section in order to show that the lower 
bound corresponding to the restricted space is sharp in the thermodynamic limit. 
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The mean-field nature of our Hamiltonian allows us to restrict the variational problem 
to a two-degrees of freedom product measures represented through the non-interacting 
Hamiltonian: 

i=l i=l 

and so, given a Hamiltonian H, we define the ansatz Gibbs state corresponding to it as 
/(cr) as: 

In order to facilitate our task, we shall express the variational principle of [H] in the 
following simple form: 

Proposition 5.3 Let a Hamiltonian H , and its associated partition function Z = 

(T 

be given. Consider an arbitrary trial Hamiltonian H and its associated partition function 
Z . The following inequality holds: 

InZ ^\nZ -0{H)+Cu{H) . (5.20) 



Given a Hamiltonian as defined in (2.2) and its associated pressure per particle pn 



^InZ, the following inequality follows from (5.20); 



N 



where 

1 



liminfpAT ^ sup Plow (5-21) 



+ahifii + (1 - a)h2fJ,2 + 

+(l-a)(-ii^ln(lii!5)-i^ln(i^)), (5.22) 



and (/ii,/i2) G [-1, 1]^ 



Proof: The (5.20) follows straightforwardly from Jensen's inequality: 
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It is convenient to express the Hamiltonian using the simbol ^ for the spins of type A 
and 7] for those of type B as: 

i,j i,j i,j i i 

indeed its expectation on the trial state is 

i,j i,j i,j i i 

(5.25) 

and a standard computation for the moments leads to 

= - a/Ar)(tanhri)2 + J^^a/N + J22((l - af - (1 - a)/A^)(tanhr2)^ + 

+722(1 — a)/N + 2Ji2a(l — 0) tanhri tanhr2) — Nahi tanhri — N{1 — a)h2 tanhr2. 

(5.26) 

Analogously, the Gibbs state of H is: 

oj{H) = —Nari tanhri — A^(l — a)r2tanhr2, 

and the non interacting partition function is: 

= ^e-^^'^) = 2^Hcoshri)^^ +2^^(coshr2)^^ 
(J 

which implies that the non-interacting pressure gives 

PAf = In Zn = In 2 + a In coshri + (1 — a) In coshr2 



So we can finally apply Proposition (5.20) in order to find a lower bound for the 
pressure = — In Z^v: 

PN = ^\nZr,^^ (^\nZr,-CuiH)+ u{H)) (5.27) 
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which exphcitly reads: 

Pat = — In Zn ^ In 2 + a In coshri + (1 — a) In coshr2 + 

+ -( Jiia^(tanhri)^ + J22(l — a)^(tanhr2)^ + 2Ji2a(l — a) tanhri tanhr2) + 
+ahi tanhri + (1 — a)h2 tanhr2 — ari tanhri — (1 — a;)r2 tanh r2 
+Jiia/2A^ + J22(l - a)/2N - Jna{tanhriY /N - J22(l - a)(tanhr2)VA^- 

(5.28) 

Taking the hm inf over N and the supremum in ri and r2 of the left hand side we get 



the (5.21) after performing the change of variables /xi = tanhri and fi2 = tanhr2. 

□ 

5.3 Exact solution of the model 

Though the functions Plow and pup are different, it is easily checked that they share the 
same local suprema. Indeed, if we differentiate both functions with respect to parameters 
/ii and we see that the extremality conditions are given in both cases by the Mean 
Field Equations: 

j2i = tanh(Jiia;yUi + Ji2(l - tt)/i2 + h) 

(5.29) 

/i2 = tanh(Ji2a/ii + J22(l - «)/i2 + h2) 
If we now use these equations to express tanh~^ yUj as a function of yUj and we substitute 
back into pup and Plow we get the same function: 

p{fii, fi2) = -^( JiiaV?+2^i2tt(l-a)AHAi2+-^22(l-tt) V2)H — ^" In -^{1-a) In 

(5.30) 

Since this function returns the value of the pressure when the couple (/xi, ^2) corre- 
sponds to an extremum, and this is the same both for Plow and pup, we have proved the 
following: 



Theorem 1 Given a hamiltonian as defined in (5.24), ^'^^ defining the pressure per 
particle as pn = — InZ, given parameters Jn, J12, J22, ^1, ^2 o-nd a, the thermodynamic 



16 



limit 

lim Pn^P 

N-*oo 

of the pressure exists, and can be expressed in one of the following equivalent forms: 

a) p= sup PLOwi/J'l, IJ'2) 

b) p^ sup PupifJ'l, IJ'2) 

6 Preliminary analytic result 

Though analysis cannot solve our problem exactly, it can tell us a what to expect when 
we solve it numerically. In particular, in this section we shall prove that, for any choice 
of the parameters, the total number of local maxima for the function p{iii, 112) is less or 
equal to five. 

We recall that the mean field equations for our two-population model are: 

Hi = tanh( Jiiq;//i + Ji2(l - q;)//2 + ^1) 
IJ,2 = tanh( Ji2Q;//i + J22(l - Q;)/X2 + /i2) 

and correspond to the stationarity conditions of p(/ii,/X2)- So, a subset of solutions to 
this system of equations are local maxima, and some among them correspond to the 
thermodynamic equilibrium. 

These equations give a two-dimensional generalization of the Curie- Weiss mean field 
equation. Solutions of the classic Curie- Weiss model can be analysed by elementary 
geometry: in our case, however, the geometry is that of 2 dimensional maps, and it pays 
to recall that Henon's map, a simingly harmless 2 dimensional diffeomorhism of M^, is 
known to exhibit full-fledged chaos. Therefore, the parametric dependence of solutions, 
and in particular the number of solutions corresponding to local maxima of p(//i,//2), is 
in no way apparent from the equations themselves. 

We can, nevertheless, recover some geometric features from the analogy with one- 
dimensional picture. For the classic Curie- Weiss equation, continuity and the Intermediate 
Value Theorem from elementary calculus assure the existence of at least one solution. In 
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higher dimensions we can resort to the analogous result, Brouwer's Fixed Point Theorem, 
which states that any continuous map on a topological closed ball has at least one fixed 
point. This theorem, applied to the smooth map R on the square [—1, 1]^, given by 

Ai2) = tanh(Jiia/ii + Ji2(l - a;)/i2 + /ii) 
-R2(/ii, yU2) = tanh(Ji2a/ii + J22(l - a)/i2 + /i2) 

establishes the existence of at least one point of thermodynamic equilibrium. 

We can gain further information by considering the precise form of the equations: by 
inverting the hyperbolic tangent in the first equation, we can /i2 as a function of /ii, and 
vice-versa for the second equation. Therefore, when J12 7^ we can rewrite the equations 
in the following fashion: 

{/i2 = , -(tanh^^ fii - Jnafii - hi) 
^12(1-") (6.31) 
Hi = - — (tanh^^ H2 - </22(l - Oi)n2 - /i2) 

Consider, for example, the first equation: this defines a function yU2(/ii), and we shall 
call its graph curve 71. Let's consider the second derivative of this function: 

5^2 _ _ 1 2/ii 
Wl ~ 'Ml -a) ' (l-/x2)2- 

We see immediately that this second derivative is strictly increasing, and that it 
changes sign exactly at zero. This implies that 71 can be divided into three monotonic 
pieces, each having strictly positive third derivative as a function of fii. The same thing 
holds for the second equation, which defines a function ^1(^2), and a corresponding curve 
72. An analytical argument easily establishes that there exist at most 9 crossing points 
of 7i and 72 (for convenience we shall label the three monotonic pieces of 71 as /, // and 
III, from left to right): since 72, too, has a strictly positive third derivative, it follows 
that it intersects each of the three monotonic pieces of 71 at most three times, and this 
leaves the number of intersections between 71 and 72 bounded above by 9 (see an example 
of this in Figure [T]) . 

By definition of the mean field equations, the stationary points of the pressure corre- 
spond to crossing points of 71 and 72. Furthermore, common sense tells us that not all of 
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these stationary points can be local maxima. This is indeed true, and it is proved by the 
following: 



Proposition 6.4 The function p{fii, fi2) admits at most 5 maxima. 



To prove 6.4 we shall need the following: 



Lemma 6.5 Say Pi and P2 are two crossing points linked by a monotonic piece of one 
of the two functions considered above. Then at most one of them is a local maximum of 
the pressure p{fii, ^12)- 



Proof of Lemma 6^: The proof consists of a simple observation about the meaning of 

our curves. The mean field equations as stationarity conditions for the pressure, so each 

of 7i and 72 are made of points where one of the two components of the gradient of 

p(yUi, ^2) vanishes. Without loss of generality assume that Pi is a maximum, and that the 

dp 

component that vanishes on the piece of curve that links Pi to P2 is — — . 

y/ii 

Since Pi is a local maximum, p{fii,fi2) locally increases on the piece of curve 7. On 
the other hand, the directional derivative of p{fii, fi2) along 7 is given by 

t ■ Vp 

where t is the unit tangent to 7. Now we just need to notice that by assumptions for any 
point in 7 t lies in the same quadrant, while Vp is vertical with a definite verse. This 
implies that the scalar product giving directional derivative is strictly non-negative over 
all 7, which prevents P2 form being a maximum. □ 



Proof of Proposition 6.4.'- The proof considers two separate cases: 



a) All crossing points can be joined in a chain by using monotonic pieces of curve such 
as the one defined in the lemma; 

b) At least one crossing point is linked to the others only by non-monotonic pieces of 
curve. 
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In 




gamma 1 
gamma 2 



Figure 1: The crossing points correspond to solutions of the mean field equations 

In case a), all stationary can be joined in chain in which no two local maxima can be 
nearest neighbours, by the lemma. Since there are at most 9 stationary points, there can 
be at most 5 local maxima. 

For case b) assume that there is a point, call it P, which is not linked to any other point 
by a monotonic piece of curve. Without loss of generality, say that P lies on I (which, we 
recall, is defined as the leftmost monotonic piece of 71). By assumption, I cannot contain 
other crossing points apart from P, for otherwise P would be monotonically linked to at 
least one of them, contradicting the assumption. On the other hand, each of // and /// 



contain at most 3 stationary points, and, by Lemma 6.5, at most 2 of these are maxima. 



So we have at most 2 maxima on each of // and III, and and at most 1 maximum on /, 
which leaves the total bounded above by 5. The cases in which P lies on II, or on III, 
are proved analogously, giving the result. 

□ 



20 



7 Comments 



The considered model generalises models which arise naturally as approximations of var- 
ious problems in theoretical physics. Furthermore, the upcoming study of social phenom- 
ena by statistical mechanics methods provides another importance source of interest for 
a model describing long-range interactions between two homogeneous populations. 

In [9] we show that it is possible to give a cultural-contact interpretation to the model 
presented here, and thanks to the mathematical results just derived, to provide non-trivial 
information about its regimes. 

It is not known at present which is the exact mathematical structure underlying social 
networks. However, it is well-accepted that interactions must be of the "small world" type 
predicted in [12], at least to some degree. We plan to return on those topics in future 
works. 

Acknowledgments. We thank Cristian Giardina and Christof Kiilske for many in- 
teresting discussions. 
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